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We examine a supersymmetric explanation for the anomalously high forward backward asymme- 
try in top pair production measured by CDF and DO. We suppose that it is due to the t— channel 
exchange of a right-handed sbottom which couples to da and ta, as is present in the R— parity vio- 
lating minimal supersymmetric standard model. We show that all Tevatron and LHC experiments' 
ti constraints may be respected for a sbottom mass between 300 and 1200 GeV, and a large Yukawa 
coupling > 2.2, yielding Afb up to 0.18. The non Standard Model contribution to the LHC charge 
asymmetry parameter is AAq — 0.017 — 0.045, small enough to be consistent with current mea- 
surements but non-zero and positive, allowing for LHC confirmation in the future within 20fb~ 1 . 
A small additional contribution to the LHC ti production cross-section is also predicted, allowing 
a further test. We estimate that 10 fb" 1 of LHC luminosity would be sufficient to rule out the 
proposal to 95% confidence level, if the measurements of the ti cross-section turn out to be centred 
on the Standard Model prediction. 
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PACS numbers: 

In the 1.96 TeV centre of mass energy pp collisions 
at the Tevatron collider, measurements of the tt forward 
backward asymmetry Afb were made by the CDF and 
DO experiments. A positive non-zero value indicates that, 
in tt production events, a higher number of events 1 N(c > 
0) had a t travelling in a more forward direction than the 
i compared to the number N(c < 0) of those events where 
the t travelled in a more backward direction than the i 



A 



FB 



N(c > 0) - N(c < 0) 
N(c > 0) + N{c < 0) 



(1) 



The Standard Model (SM) prediction of this quantity is 
A F *£ = 0.066 ± 0.020 and dominantly derives from 
the interference between tree-level and one-loop quantum 
chromodynamics (QCD) diagrams. CDF measured an 
unfolded value 2 0.158±0.075 [3,11], whereas DO measured 
0.196±0.065 155, each significantly higher than Apg. We 
suppose here that the discrepancy between measurements 
and the SM prediction for Afb is due to a particular 
beyond the SM process, and examine other constraints 
to see if the explanation remains viable. 

Since the LHC is a pp collider and thus has an ini- 
tial state which is symmetric under c O — c, A F ^ C = 0. 
However, the LHC is able to measure a related but dif- 
ferent charge asymmetry in the number of tops that are 



c = cos 8, where 9 is defined to be the scattering angle between 
incoming proton beam and outgoing top in the centre of mass 
frame of ti. 

Different shower models can produce different values when they 
are used in the unfolding because of the different treatments 
of QCD coherence Q, yielding effects on Apg of order several 
percent. 



travelling at a smaller angle to the beam-line compared 
to the number of anti-tops that are travelling closer to 
the beam-line: 
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where = l/21n(£' i + p iz )/{Ei — p iz ) is the rapid- 
ity of particle i. The SM prediction for 7 TeV col- 
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0.006 ± 0.002, which is consistent 
with the combined ATLAS and CMS measurements 
A v c = -0.015 ± 0.04 @, 0. Any non-SM explanation 
for the high measured value of Apg must also therefore 
not predict too high a value for a non-SM contribution 
AA V C = A V C - A y c SM . 

Many models have so far been proposed to explain 
Apb measurements and several have failed other con- 
straints. Axigluons, W and Z' vector bosons, as well as 
the t— channel exchange of various scalars [8H10l| have all 
been proposed. Here, we show that the exchange of a par- 
ticular scalar: a right-handed sbottom, can explain the 
apparent sizable enhancement to Afb while respecting 
other relevant constraints. Ref. [ll| included this, among 
other possibilities, as an explanation for Afb- However, 
the authors only considered couplings smaller than 1.25 
because they required perturbativity up until the GUT 
scale, and found that the new physics contributions to 
Afb were too small to obtain it to within the la mea- 
surement errors, although they were able to obtain it to 
just within 2a measurement errors. We shall consider 
larger values of the coupling, which we show are neces- 
sary to explain the data. Our work also goes further than 
Ref. |ll| in the sense that we consider the LHC charge 



asymmetry and top production constraints as well. 
The R— parity violating (RPV) interactions of the min- 
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imal supersymmetric standard model (MSSM) include 
the superpotential term 



W = 



313 



TrDrBr, 



(3) 



where gauge indices have been suppressed and Tr, Dr, 
Br are chiral superfields containing the anti-right handed 
top tR, anti-right handed down d,R and anti-right handed 
bottom quark 6r, respectively. W — \'^ 12 TrDrSr could 
also explain Apr, although the constraints on rris R are 
likely to be stronger than those on the Sr, since the 
strange PDFs are higher than the bottom PDFs, and 
so they will be predicted to be produced more readily at 
the LHC. The operator in Eq. [3] has an antisymmetric 
colour structure. It leads to an additional tree-level pro- 
cess that contributes to Apr, the Feynman diagram of 
which is shown in Fig. [T] Ref. llj put an upper bound 



on A 313 < 1.25 on the grounds of perturbativity up to 
the GUT scale. In the present paper, we shall not worry 
about a premature ultra-violet completion of our model: 
the coupling will reach a Landau pole around 10-100 TeV 
or so, and above that scale the effective theory will then 
change. We shall consider couplings as large as 5, al- 
though we note that at the higher end above 3.5, our pre- 
dictions may start to become less accurate due to larger 
higher loop corrections that we do not take into account. 

The tree-level A 313 ^ contribution to the differ- 
ential cross-section of tt production is, where T = 



IX" I 2 
l A 313 



dAa 



m 2 /m 2 ) 2 /(8ir) is the sbottom width, 
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where /3 = \f\ — 4mf / s, s = (jp\ +P2) 2 , ot s is the strong 
coupling constant and rrit is the top quark mass. Eq. [3] 
predicts a non-zero contribution to Apr since it is not 
even in c. It disagrees with the expression for dAa /dt in 
Ref. 12[ by a factor of 64 in the first term, but agrees with 
the recent Ref. [Hf]. In order to calculate the non-SM 
contribution to the asymmetry A Apr, we must convo- 
lute the differential cross-section with the parton distri- 
bution functions (PDFs) for (anti-)down quarks in (anti- 
protons numerically. We have numerically checked Eq.0] 
at a test point for large A 313 (where the first term dom- 
inates) that integrating Eq. [4] over the PDFs gives an 
identical result to the program MadGraph.5_vl_4_5 [14\ . 

Of eight scalar models considered in Ref. that 
might have explained the Apr measurements, the ex- 
change of a charge —1/3 colour triplet scalar that has a 
large coupling to dRtR was considered and discarded. As 
far as Apr goes, this is identical to our SUSY model ex- 
cept for the other interactions of the MSSM (which are 



important for passing the atomic parity violation con- 
straints, as mentioned below). It was deduced that for 
this scalar, there is no parameter space that simultane- 
ously satisfies other constraints as well as AA FB > 0.2, 
where implies the high invariant mass constraint m t j > 
450 GeV. However, CDF data on AA FB has lowered 
considerably since Ref. (loj . with a new CDF measure- 
ment implying AA h FB = 0.20 ± 0.07. We shall show 
that the model now agrees with all current relevant mea- 
surements. The measured value A l FB = — 0.116 ± 0.153, 
which is the asymmetry with the additional cut m t i < 
450 GeV, is consistent with the Standard Model predic- 
tion A l FB SM = 0.04. 

We shall work in a 'bottom-up' framework, where 
we deal with weak-scale SUSY, making no assumptions 
about the ultra-violet limit of our model (in particular, 
we do not demand perturbativity of couplings up to a 
grand unified scale). Thus, for the present, we shall work 
in a simplified model where we set all sparticles of the 
MSSM to be heavy, aside from the right-handed sbottom 
bR, unless for some particular reason we require another 
sparticle to play a role. The coupling A 313 , if it is the only 
non-negligible real RPV coupling, is not constrained to 
be small by indirect experimental data [HI, [13] for 1>rs 
that are not too light 3 . The stron gest constraints from 
Ri = T(Z a hadrons)/r(Z° 11) [HI are still too weak 
(for the right-handed sbottom masses > 300 GeV that 
we shall be interested in) to be limiting. Eq. [3] has ex- 
actly the right properties to evade current stringent flavor 
constraints on di-quarks: it only couples to right-handed 
quarks and induces no tree-level flavour changing neu- 
tral currents [2(|. Entertaining the possibility of other 
baryon-number violating couplings X'[^ k in the presence of 
a large coupling A 313 , there is a particularly strict bound 
on |A 3 ' 13 A 3 ' 23 *| < 0.01 [ll| coming from K° - K° mixing 
constraints for sparticle mass input parameters less than 
1 TeV. Thus, strong constraints upon |A 3 ' 23 | would ap- 
ply, and some fermion mass model building would have 
to be performed to see if this were realistic. It has been 
argued that simple Higgsed abelian flavour symmetries 
and the Froggatt-Nielsen mechanism would naively pre- 
dict that the order of magnitude of |A 3 ' 23 | would be larger 
than that of I A 313 1 2l|. However, this order of magnitude 
estimate is rather rough, and assumes that several cou- 
plings involved in the Froggatt-Nielsen mechanism are 
of order 1 , an assumption that may be violated in more 
explicit constructions. Here, it seems premature to ap- 
ply any such model building, since as we shall show, LHC 



3 If we were to have light charginos and left-handed down squarks, 
as well as a mixing between left and right-handed down squarks, 
we would have severe constraints upon Ag 13 coming from a loop 
diagram inducing neutrino anti- neutron oscillations Hal . For in- 
stance, if all sparticles were to have a mass of 200 GeV (and the 
left and right-handed squarks were to be mixed with a trilinear 
scalar coupling of 200 GeV), then the upper bound would be 
A" < 0(0.04). 
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FIG. 1: Leading order supersymmetric (SUSY) contribution 
to Afb from the RPV coupling A313 and right-handed sbot- 
tom exchange. Pi,P2,Qi and Q2 denote 4- momenta. 



0.037 < AA FB < 0.205 
-0.65 < Aa% BV /ph < 1.51 



-19.2 < Aatf C7 /pb < 39.2 



-0.079 < AA y c < 0.061 
-76 < A<T^- £;v '(bin)/fb< 76 
0.062 < AAp B < 0.33 



TABLE I: 95% CL constraints on new physics contributions 
to observables that are brought to bear upon our model. The 
limits have been derived by using naive summation in quadra- 
ture of all errors. 



charge asymmetry and top production cross-section mea- 
surements can test our proposal. We thus proceed on the 
basis that |A 323 | is small enough to evade this bound, ap- 
proximating it by zero in our numerical analysis. 

We now discuss the various constraints that we shall 
employ. Naively adding errors in quadrature, we find a 
CDF and DO weighted average of A FB = 0.187 ± 0.037. 
When combined with the SM prediction, this implies 
a measured non-SM contribution of A Afb — Afb — 
A S F M = 0.121 ±0.042. On the other hand, Tevatron mea- 
surements of the SM total tt production cross-section are 
roughly in line with SM predictions. A non-SM com- 
ponent is restricted by 22] AaT f EV = 0.43 ± 0.54 pb. 



Aojfv 

The differential production cross-section da^ EV jdr 



was measured by CDF in Ref. 



We shall employ 



(following Ref. [10|) the measurement o^ w (700 GeV < 
m t{ < 800 GeV) = 80 ± 37 fb, versus a SM predic- 
tion of [H, Hi| 80±8 pb, so the non-SM contribution 
A<r£ (bin) must not be too large. This invariant mass 
bin is far away from the bulk of the tt differential cross- 
section, and so ought to provide information that is ap- 
proximately independent of the information from crJi EV ■ 
At the 95% confidence level (CL), this leaves little room 
for a non-SM contribution of Aaf t EV (bin) . ATLAS and 
CMS have [H, H3 measured the 7 TeV pp tt cross- 
section to be cr^ Bcr = 173.4 ± 10.6 pb, versus a SM 
prediction of 163 ± 10 pb. In Table [fl we summarise the 
constraints that we require predicted observables to sat- 
isfy. 

Throughout the present paper, we calculate experi- 
mental observables using the matrix element event gener- 
ator MadGraph5_vl_l_4_5 [II] assuming m t = 173.1 GeV, 
the CTEQ6L1 PDFs [H| and using the FeynRules (29| im- 
plementation of the RPV MSSM [U |3l|. We define an 



11 by 11 grid in -X'^ 13 parameter space, simulating 
100000 tt production events and interpolating predicted 
observables in between the grid points. 

Fig. [5] shows the effect of the constraints upon the pa- 
rameter space of the model. Encouragingly, there is some 
parameter space which simultaneously fits all of the rele- 
vant current tt constraints, shown by the region enclosed 
by the solid contour. By drawing additional contours 
which we do not show here for reasons of clarity, we find 
the values of the various observables which are within 
the good fit region. These predictions are displayed in 
Table [TT] and should aid future tests of the model. AT- 
LAS obtained a statistical error on the charge asymme- 
try of 0.028, and a systematic error of 0.024 in 1 fb— 1 
of 7 TeV LHC data We expect both errors to scale 
like 1/Vc, since the systematics are also measurement 
dominated. Thus we expect, for 10 _1 fb, that each error 
will be around 0.01, which will certainly test most of the 
parameter space, where AA V C is higher. 

On the day of completion of this article, Ref. [HJ ap- 
peared on the electronic archive. It considers the effect 
of the exchange of a charge —1/3 scalar triplet diquark, 
and the bounds on parameter space coming from Afb, 
Ac, vJi^ and o-^ HC7 measurements agree with those 
shown here and aside from the fact that the paper also 
uses MadGraph to do the simulations, it provides an in- 
dependent confirmation of some of our results. 

All of the predictions above apply to a generic model 
containing a colour triplet scalar of charge —1/3 coupling 
to dtf tfi . Such models are extremely strongly constrained 
by atomic parity violation (APV) constraints [33]. The 
RPV MSSM, however, has many additional particles and 
interactions and the possibility for cancellations of differ- 
ent contributions to the APV is open. Indeed, a very 
recent paper (l3l | has shown that stop mixing contri- 
butions to APV can completely cancel the contribution 
from A 313 ^ 0. 

It was recently argued that baryon number violating 
couplings such as A 313 and a light stop could simultane- 
ously explain the naturalness of SUSY and its evasion 
of 7 TeV LHC searches based on large missing trans- 
verse momentum [34] ]. SUSY implies additional inter- 
actions with identical coupling strengths: in particular, 
£ = (A 3 i 3 )*tfl'iflC' t 6ij + B.C. + . . ., where C is the charge 
conjugation matrix and T denotes transpose. If tfj is not 
too heavy, t— channel stop exchange should induce a bb 
forward backward asymmetry [35| . One may also expect 
resonant stop production to produce a bump in the rrijj 
distribution in two jet events where one of the jets is a 
b. However, such a signal will be suppressed by small 
bottom PDFs and further study is required to establish 
the viability of detection. 

In summary, we have shown that the RPV MSSM 
can explain the anomalously large Ap b measured at the 
Tevatron experiments, provided that couplings are cho- 
sen in a region that may violate perturbativity below the 
GUT scale. The t— channel exchange of a right-handed 
sbottom of mass 300-1200 GeV which couples to dB.tR 
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FIG. 2: Predicted non-SM contributions to various constraining observables: (a) the non-SM ti forward backward asymmetry 
parameter, (b) the total ti Tevatron cross-section, (c) the total ti Tevatron cross-section in the bin 700 GeV< m t i <800 GeV, 
(d) the LHC charge asymmetry. Each individual constraint is respected below the broken contour, except for in (a), where it 
is respected above the contour. Inside the solid contour, all limits in Table U are respected. We have checked that the 8 TeV 
LHC has the same predictions for A v c as those displayed in (d), to a very good approximation. 



0.037 < AA FB < 0.11 
-0.1 < Aa^ EV /pb < l.i 
0.006 < AA l FB < 0.07 
8 < Aa^ HC7 /pb < 25 



0.017 < AA y c < 0.045 
!5 < Aa£- BV (bin)/fb< 76 
0.062 < AA h FB < 0.16 
13 < Aa^ HC8 /pb < 33 



TABLE II: Predicted values of new physics contribution to 
various observables for parameter space that passes all current 
experimental constraints. 



with an interaction strength bigger than 2.2 is required. 
There is parameter space which passes all relevant con- 
straints from the Tevatron and the LHC. Future LHC 
signals are predicted to be: a small positive asymmetry 
parameter 0.017 < AA V C < 0.045 and a ti production 
cross-section that should be larger than the SM predic- 
tion by at least 8 (13) pb at 7 (8) TeV. Assuming that 



the uncertainties oc 1/ v C (where C is the integrated lu- 
minosity), there should be enough information in 10 fb -1 
of 7 TeV LHC data to exclude such an enhancement of 
the production cross-section to the 95% confidence level. 
Further confirmation through the top charge asymme- 
try may take longer, but is within reach of 20 fb . 
These predictions, along with the other predictions in 
Table UH provide targets for ongoing LHC tests of the 
model. Top spin observables also provide an interesting 
way to discriminate models of new physics that explain 



the Tevatron measurements Ap b [36l . l37j , and it would 



be an interesting future project to investigate the predic- 
tions for them coming from X'^ 13 . If the sbottoms are not 
too heavy, then there is the possibility of confirming the 
mechanism through bab* R production at the LHC, which 
would then decay to ti plus 2 jets, with an invariant mass 
bump in the combined t and jet mass. A discovery such 



5 



as this would provide a definitive test of the model. We 
leave its investigation to future work. 
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